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We calculate the quasiparticle contribution to the zero temperature Hall conduc- 
tance of two-dimensional extreme type-II superconductors in a high magnetic field, 
using the Landau basis. As one enters the superconducting phase the Hall conduc- 
tance is renormalized to smaller values, with respect to the normal state result, until 
a quantum level-crossing transition is reached. At high values of the order parame- 
ter, where the quasiparticles are bound to the vortex cores, the Hall conductance is 
expected to tend to zero due to a theorem of Thouless. 



I. INTRODUCTION 

Numerous superconductors like high-T c cuprates, Al5's, boro-carbides and many organics 
reveal the importance of Landau quantization at high magnetic fields and low temperatures, 
as evidenced, for instance, by the recent observation of the de Haas-van Aiphen oscillations 
in the superconducting phase |]. In these systems the cyclotron splitting (u c ) between 
the Landau levels (LL) is the largest energy scale, in particular u c > A (where A is the 
superconductor order parameter). The effect of the Landau quantization has attracted 
theoretical attention and reveals itself by the appearance of gapless superconductivity |]2]|| 
and power law dependence of the low-T thermodynamic properties. A good description 
of the de Haas-van Aiphen oscillations has been obtained taking into account the Landau 
quantization [|J. The gapless nature of the spectrum also leads to an enhanced accoustic 
attenuation in contrast to usual gapped superconductors H and it survives even in the 
presence of disorder |J. 

Close to the upper critical field, and at small values of A, a diagonal approximation 
(DA), where the coupling between the LL is neglected, is appropriate [^,0. In this regime 
the quasiparticles propagate coherently throughout the vortex lattice. Associated with the 
zeros of the order parameter in real space are gapless points in the magnetic Brillouin zone. 
As A grows, the coupling between the LL has to be taken into account. A perturbation 
scheme in the off-diagonal couplings can then be carried out, as long as there are no level 
crossings 0, and it can be shown analytically to all orders in the perturbation theory that 
there is always a discrete set of gapless points. In general, a numerical solution of the 
Bogoliubov-de Gennes (BdG) equations has to be carried out. For large values of A a 
regime is reached where the quasiparticles are bound to the vortex cores with an energy 
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spacing typical of an isolated single vortex ||. In this regime one expects that a tight- 
binding approximation should yield good results, since in this low-field regime the vortex 
cores are sparse. Such a procedure has been carried out and a gapped regime is indeed found 
IIH signalling a transition from the high-field gapless regime to the low-field gapped regime. 
Between the small A regime and this gapped regime several quantum level-crossings take 
place. This is a difficult numerical problem because the number of Landau levels increases 
considerably as the field lowers and it is not easy to determine the value of A for which the 
transition to the gapped regime takes place. 

On the other hand, Thouless has shown that the Landau levels do not constitute a 
complete basis to be used in a tight-binding approximation scheme and therefore it is not 
possible to construct well localized representations of the magnetic translation group unless 
the Hall conductance, a xy , is zero [TO]. Thouless showed that if Wannier functions are 



constructed they decay with distance at most as \r — Ri\~ 2 (if the Hall conductance is 
nonvanishing) but he also showed that well localized Wannier functions (decaying faster 
than any power) can be constructed if a xy vanishes. Using the Balian-Low theorem it has 



also been shown that in general | L] the Landau levels cannot be chosen sufficiently localized 
to make the Ax and Ay uncertainties finite. This is equivalent to the statement by Thouless 
on the slow decay rate of the Wannier functions with distance. The notable exception is 
a sub-band with zero Hall conductance. Since it is intuitive that in the limit the vortex 
concentration is low a tight-binding description should be appropriate, we expect, in the 
light of Thouless theorem, that a xy should be zero in this low-field regime. This is suggested 
by the localized nature evidenced by the numerical solution of the BdG equations. To 
determine if the Wannier states are indeed a good basis is a delicate and difficult matter and 
in general a tight-binding description will have one conducting state at the band center [FTS 



We expect however that in this low-field regime the functions will be well localized and we 
argue therefore that the Hall conductance could be used as an order parameter to signal the 
transition from the high-field gapless regime, where Landau quantization has been shown to 
occur, (finite a xy ) to the low-field gapped regime (zero a xy ). 

The calculation of the Hall conductance in the superconducting phase is also interesting 
by itself. Several authors have paid attention to the Hall conductance of the vortex lattice. 
In general, there are two contributions: one is due to the vortex motion and the other to 
the quasiparticle contribution (usually associated with modes localized in the normal region 
inside the vortex cores). One of the reasons for this interest is that in the superconducting 
phase a xy has a different sign with respect to the one of the normal phase ||13|| . This has 



been shown to be due to the vortex motion part since the localized modes are predicted to 
give a contribution with the same sign as the normal phase value. 

In this paper we focus solely on the quasiparticle contribution which is the only remaining 
if the vortex lattice is pinned to some imperfection (in this case the vortex motion is frozen). 
We do not address here the sign change but study the influence of the coherent propagation 
of the quasiparticles. Also, we take into consideration the Landau quantization in a regime 



where a quasiclassical approximation is not valid [14 



II. CALCULATION OF THE HALL CONDUCTANCE 

We calculate the Hall conductance using the Kubo formula 
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a xy (r,r') = -ihL 2 £{< 0\J x (r)\f3 >< (3\J y (r')\0 > 

- < 0\J y (r')\f3 >< {3\J x (r)\0 >} 1 (1) 
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where the currents are given by 
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Here a is the spin projection, i — x, y, z, A is the vector potential, Vv( r ) is the electron field 
operator, the energies ep are the full solution of the many-body problem and (3 = is the 
groundstate. We consider a square lattice of side L. 

We calculate the Hall conductance obtaining the energies from the solution of the BdG 
equations [0,13]]. The field operators are written as 



^W = EK( r )7,, q , T -^(r)7:, qa ; 
^W = E(<( r )7,, q a + <(r)7l, q , T ) (3) 

and 

6 /3 = E e q7i, q ,a>,q, CT - (4) 



tr,z/,q 



Here 7j qcr creates a quasi-particle in the level u, with momentum q and spin a. The 
problem is diagonal in q (vectors of the magnetic Brillouin zone) due to the periodicity of 
the Abrikosov vortex lattice. The amplitudes u and v are expanded in the Landau basis (n 
is the Landau level) as 

<( r ) = XXrAn(r) (5) 



and 



<( r ) = EC^ q nW. (6) 



Here qn (r) are the eigenfunctions of the magnetic translation group in the Landau gauge 
(A x = —Hy, A y = 0, A z = 0) belonging to the n th Landau level. The amplitudes u^ n and 
i>q n are the solutions of the BdG equations and are the energy eigenvalues. The BdG 
equations to be solved are 
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where e n = (n + 1/2) — h/Tioj c (/i is the chemical potential), A nm (q) is the matrix element 
of the order parameter A(r) between electronic states (q, n) and (— q, m) 0]. In general, 
the order parameter is expanded in the LL for a charge 2e, which leads to a complicated 
structure for the BdG equations which have to be solved self-consistently, since the order 
parameter satisfies A(r) = V < V ; t( r )V ; j.( r ) >j where V is the electron-electron interaction. 
In the DA these equations are easilly solved 
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(and therefore, u = n). In this case 

<n«-q» = (9) 
qn qn o,n v ' 

As A grows, the off-diagonal terms become increasingly important. Recently, the leading 
corrections to the normal and pairing self-energies have been found by delaying the Nambu 
rotation to the last step and it was obtained that the renormalized quasiparticle energies 
can be written as in the DA like 



e^ q = ±Ve n (q) + |A nn (q)| 2 (10) 
but with renormalized normal and pairing terms which to leading order are given by 

e„(q) = e n + ^ 

A nn (q) = A„ n (q) 

_ y-^ ^ n > n +P (O) ^n+p,n+p' ( ~ C l) ^n+p' ,n (q) 
p^0,p'^0 PP 

The effect of the off-diagonal couplings is taken into account by renormalizing the diagonal 
(in the Landau index) terms keeping a diagonal problem to be Nambu rotated (y is still 
v = n). 

In Fig. la we show the energy of the lowest band above the Fermi level, obtained from 
eq. (8), as a function of q (we only show one quadrant of the magnetic Brillouin zone due 
to the square lattice symmetry) for A = 0.1 and n c = 10 (we fix \i at this LL energy). In 
Figs, lb and lc we show the same band taking into account the off-diagonal terms in eq. 
(7) considering A = 0.1 and A = 0.5, respectively. Results similar to these were previously 
presented in other forms (see e.g. @,§) but we present them here as well to stress the role 
of the off-diagonal terms. The energy of the lowest band above the Fermi level is given by 
the gap function |A nn (g)| in the DA. This is given by || 
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A ™(D = J2e 2ik ^-^ k ^ l2 H 2n [V2(q x + t3)Z] (12) 



nl 



where a is the lattice constant of the square vortex lattice and / is the magnetic length, 
/ = \J~^ = -7=. The gap |A nn (g*)| has zeros in the magnetic Brillouin zone in momentum 
points (qj) which are in direct correspondence with the positions of the vortices in real space 
(z^ such that qjl = Zi/l 0. These gapless points are due to the center of mass motion of the 
Cooper pairs in the presence of a high magnetic field and not due to the internal structure 
of the gap function as in (i-wave superconductors. The gapless points grow considerably in 
number as n c grows 0. However, they preserve in the magnetic cell the vorticity per unit cell 
of the order parameter in real space. Away from the gapless points the gap raises to a scale 
of the order of A giving rise to a complicated structure that oscillates ever more strongly 
as n grows. Increasing A only rescales the energy. The energy has a linear dispersion 
around most of the zeros (—1 vorticity). Other sets of zeros have a quadratic dispersion 
for the square (triangular) lattice corresponding to a —2 (—3) vorticity, respectively. This 
dispersion leads to power law behavior at low temperatures, as mentioned above. Including 
off-diagonal terms small gaps open in the spectrum |7j which are due to the normal (pairing) 
part of the self-energy and which are of order A 2 (A 3 ). However, throughout the Brillouin 
zone the numerical prefactors are very small and the spectrum is similar almost everywhere 
to the one obtained in the DA. The notable exception are the Eilenberger points (EP) (e.g. 
q = (—0.5,0.5)) where the gap increases significantly as A grows 0. The off-diagonal 
matrix elements A nm (q) have in general a different structure of zeros. If n + m is even the 
Eilenberger points remain but otherwise they are absent. 

These results are best summarized considering the density of states (DOS). In Fig. 2 
we plot the DOS as a function of energy for different values of A. For small A the DOS 
is broadened from the Landau level locations by an amount of the order of A. The DOS 
is high at low energies. As A grows the energy interval grows and eventually at A of the 
order of the cyclotron energy the lowest band and the next band approach each other and a 
quantum level-crossing transition occurs. Note that as A grows and the DOS spreads over 
to higher energies, there remains a high DOS at low energies due to the low gap states that 
retain the characteristic spectrum of the DA. Therefore, it is expected that the DA will be 
appropriate at low A 0,[J (see however ref. || for a thorough discussion). 

To calculate the Hall conductance we insert eqs. (2,3) in the expression for a xy eq. (1) 
and obtain for the average Hall conductance 
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— J dx J dx' J dy J dy'a xy (x,y; x' ,y') (13) 



after several straightforward but lengthy steps (q = 2e) 



+ bV<(<n + l^T + <n-l\f^)( V -\rn + ij^- «^m-l^)]} (14) 
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The sum is extended over the excited states. Since the current operator is a two-particle 
operator the excited states only differ from the groundstate by the occupation of two quasi- 
particle (fermionic) states (this is the reason for the prime in the number occupation term). 
Here N<p is the number of momenta appearing in the sum which equals the number of vortices 
in the system. The energies have been rescaled by hu c . 

A selection rule in the DA implies that v 7^ v' . We can therefore sum over the spin 
variables. The expression for a then takes the form 

a 1^ 1 1 



i 4 w ^ - e o) 2 n* 



K{(n„ q n n+ i iq )'(n+ 1)[ 



A* nn A n+1[n+1 _ e n + e n+ i 

q q q c q 

,n - ^n— 1 



A* A ,n . ,n—i _i_ 

-K^J'M^^Fi - "~ " - + (15) 



q c q c q q 



Here €p = e™ + e™ +1 in the first term and €p = e™ + e™ -1 in the second term. Fixing the 
chemical potential at the level n c (n c + 1 occupied levels) and taking A = (normal phase) 
we get that a = ah/q 2 — n c + 1, as expected. The only contribution comes from the first 
term with n = n c . Taking now a small value of A the two dominant contributions come from 
the previous term and from the first term with n — n c — 1. This leads to a discontinuity in a: 
at A = the Hall conductance a = (n c + l) while at small, but nonzero A, a ~ (n c +l) — 1/2. 
As A grows a decreases continuously. 

In Fig. 3 we show a as a function of A for the DA and for the leading order perturbation 
theory (PT) in the range of values of A up to the order of the first level-crossing. We 
consider the cases n c — 4, 10. As A grows the off-diagonal terms renormalize downwards 
the Hall conductance with respect to the DA value. In the normal phase (A = 0) the Fermi 
level (for a completely filled level) is in the gap between two Landau levels. As A is turned 
on, /i is kept fixed at /i = n c + 1/2 + i] (rj — > 0) and the levels spread to higher energies which 
leads to a decrease of the Hall conductance since the fraction of low energy states that may 
conduct decreases. The presence of the off-diagonal terms fastens the rate of decrease for 
larger values of A, for the same reason. At very low A the two methods agree, as expected. 
Close to the level crossing the study is very difficult because a fully self-consistent calculation 
is needed to properly cover the crossover to the gapped regime. Going beyond the DA the 
several eigenfunctions of the BdG equations have components in the various Landau indices, 
which due to the gapless (or almost gapless) nature of the spectrum causes difficult numerical 
problems due to the energy denominator in eq. (13). Also, strong mixing of the LL destroys 
the LL structure beyond the level-crossing(s) transition(s). There is a discontinuity in the 
Fermi level due to this level crossing (by ±2u c because of the doubling due to the particle 
and hole bands, u and v, respectively). After the transition, a is expeced to decrease again 
and to tend to zero as A grows even further (eventually after several level-crossings). 



III. CONCLUSIONS 

In summary, we have calculated the Hall conductance of a pinned vortex lattice in a 
high magnetic field using the Kubo formula and the solution of the Bogoliubov-de Gennes 
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equations for the wave function amplitudes expanded in the Landau basis. We compared 
the diagonal approximation with the leading order perturbation theory recently introduced 
0. The Hall conductance decreases from the normal state value due to the presence of 
the low-lying states immediately above the Fermi energy. As A grows the spread in energy 
increases and a decreases. We limited the study to the region of validity of the PT. We 
suggest that a may be used as an order parameter to detect the transition from the gapless 
regime (finite a) to the gapped region (zero a) where a tight-binding description should be 
appropriate. 

The author acknowledges discussions with Zlatko Tesanovic at the early stages of this 
work and for calling his attention to ref. 10. This work was partially supported through 
PRAXIS Project /2/2.1/FIS/302/94. 
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FIGURE CAPTIONS 

Fig. 1- Lowest energy band as a function of q for a) A = 0.1 in the DA and including 
off-diagonal terms for b) A = 0.1 and c) A = 0.5, respectively. Only one quadrant of the 
Brillouin zone is shown. 

Fig. 2- Density of states (DOS) in arbitrary units as a function of energy for n c = 10. 
The first two bands are shown. 

Fig. 3- Hall conductance for n c = 4, 10 as a function of A for the two methods considered: 
diagonal approximation (DA) and leading order perturbation theory (PT). The two points 
at n c + 1 are the A = result (normal phase). 
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